Relativistic Tunneling Through Two Successive Barriers 
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We study the relativistic quantum mechanical problem of a Dirac particle tunneling through two 
successive electrostatic barriers. Our aim is to study the emergence of the so-called Generalized 
Hartman Effect, an effect observed in the context of nonrelativistic tunneling as well as in its 
electromagnetic counterparts, and which is often associated with the possibility of superluminal 
velocities in the tunneling process. We discuss the behavior of both the phase (or group) tunneling 
time and the dwell time, and show that in the limit of opaque barriers the relativistic theory 
also allows the emergence of the Generalized Hartman Effect. We compare our results with the 
nonrelativistic ones and discuss their interpretation. 

PACS numbers: 03.65.Xp, 03.65.Pm, 73.40.Gk 



I. INTRODUCTION 

The phenomenon of tunneling is one of the most strik- 
ing and extensively studied consequences of quantum me- 
chanics. Yet,after decades of scrutiny (for reviews see, 
for example, P, 0,11,13]) it still presents serious conceptual 
challenges, such as a meaningful definition of tunneling 
time, that is, the time it takes for a particle to tunnel 
through a potential barrier. 

Several different scales of time associated with the tun- 
neling process have been proposed (see Among the 
most prominent ones are the phase time (or group delay 
time) , given by the energy derivative of the phase shift in 
the transmission (or reflection) amplitude, and the dwell 
time, which is related to the average time spent by the 
particle in the region of the potential. 

It is well known that for the tunneling of a particle 
through an opaque barrier the group delay saturates with 
the width of the barrier, a phenomenon called Hartmann 
effect [5| . Several authors interpret this saturated time as 
the transit time for the particle to go through the poten- 
tial, which would imply, as an immediate consequence, 
the possibility of superluminal (group) velocities for bar- 
riers with a sufficiently large spatial extension. Such an 
interpretation has been in the center of an intense debate 
in the literature (see, e.g., 0, H, 0] and references there 
cited). 

Recently an apparently even more paradoxical effect, 
which became known as Generalized Hartman Effect, has 
been brought to attention, not only in the context of non- 
relativistic quantum tunneling, but also in the context of 
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its electromagnetic counterparts. This effect consists in 
the fact that for tunneling through two potential barriers 
separated by a distance I the phase time is, in the limit 
of opaque barriers, independent not only on the barrier 
widths but also on the spacing between them Q (see also 
Efl|)- In fact, Esposito [ll| showed that for a sys- 
tem of N barriers the phase time is independent also on 
the number of barriers. Despite the fact that phase time 
cannot, in general, be interpreted as a propagation (or 
transit) time for the particle (or wavepacket), this effect 
is counterintuitive since one would, naively, expect that 
in the space between the barriers the group delay could 
be viewed as a propagation time, and therefore, it should 
depend on the distance between the barriers. 

In the last years several papers have also analyzed the 
problem of q uantum tunne ling from a relativistic stand- 
point [l| El [3 [H, [H, El EI [H. Some of these papers 
were concerned with the fact that the analysis of possible 
superluminal group velocities associated with the Hart- 
mann effect should be properly addressed in the context 
of a relativistic theory [H El El El EH ; others were 
concerned with general aspects of the relativistic prob- 
lem, such as the relation between phase time and dwell 
time 11811 or the relation between dwell time and Larmor 
times In fact, due to the relevance of the tunnel- 

ing phenomenon, it is important to consider the possible 
quantitative and/or qualitative differences in the phase 
(and dwell) time arising due to the relativistic dynamic. 
What is more, a clear understanding of the relativistic 
aspects of tunneling is imperative if one wants to eventu- 
ally obtain a meaningful time scale for this phenomenon, 
because the instantaneous spread of the probability den- 
sity in nonrelativistic quantum mechanics (2p| makes it 
difficult to define an unambiguous tunneling time in the 
context of Schrodinger theory. 

Most of the above papers were concerned with a sin- 
gle potential barrier or potential well. In fact, to the 
best of our knowledge, the only relativistic treatment of 
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the two barrier case was due to Leavens and Aers [21|, 
which were concerned with Larmor times at resonance. 
Thus, the present work complements the previous ones 
by considering the relativistic approach to the problem 
of a wave packet tunneling through two successive bar- 
riers. We shall address the emergence of the generalized 
Hartman effect in this context and discuss its possible 
interpretations. We also compare our results with those 
obtained in the nonrelativistic framework. 
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II. THE RELATIVISTIC PHASE AND DWELL 
TIMES 

Here we shall be concerned with the relativistic one- 
dimensional scattering of a mass m and spin- 1/2 wave 
packet by an electrostatic time-independent potential 
V(z). The general form of the incident wave packet is 
given by 



V(z,t) = / dE A(E)ijj{z)( 



(1) 



where A(E) are its Fourier coefficients, while ip(z) must 
satisfy the time-independent Dirac equation associated 
with the energy E, 

{-ihca z d z + (3mc 2 + V{z)} ij)(z) = E^{z) , (2) 

with a z and (3 being the Dirac matrices (we will follow 
the conventions in 22]). The potential in which we are 
interested consists of two potential barriers of height Vb 
and width a, and spaced by a distance I, as seen in Figure 
[1] (since barriers of different heights and widths do not 
introduce any novelty, they will not be considered here). 

Moreover, we shall consider an incident wave packet 
whose Fourier energy distribution A{E) is very sharply 
concentrated around a given value Eq, corresponding, 
therefore, to a smooth modulation of the cigenfunction 
corresponding to Eq. Such a wave packet must have, 
therefore, a large spatial extension. For our purposes we 
shall assume that the energy dispersion of the wavepacket 
is sufficiently narrow such that its spatial extension is al- 
ways very large when compared to the extension of the 
region in which the potential is nonvanishing (the re- 
gion < z < I + 2a in Figure [T]). With these assump- 
tions in mind it is justifiable the use of the stationary 
phase method to follow the position of the peak of the 
wavepacket in the free regions (regions / and V) (see, 
for example, [23| and references therein). We shall also 
consider that Eq is a positive energy (particle) in the 
evanescent region, Eq — mc 2 < Vb < Eq + mc 2 . There- 
fore, the region of supercritical potential, in which there 
is pair production (and the associated Klein paradox) and 
where, therefore, the one-particle Dirac equation ceases 
to be valid, will not be considered (for a study of the su- 
percritical region for the one barrier potential see [2~jj). 

Considering, as usual, a wave packet incident only from 
the left and having spin up (we make this assumption 



FIG. 1: Schematic view of the two barriers arrangement. We 
consider two rectangular identical barriers of heigth Vb and 
width a, separated between themselves by a distance L. 



without loss of generality, because the potential consid- 
ered here causes no spin flip), the general solution of the 
stationary problem in the various regions indicated in 
Figure Q] is given by 

iPj(z) = e lkz u E {k) + Re- lkz u E {-k) , 

ipn[z) = Ae- qz u E - Vo (iq) + Be qz u E ^ Vo (~iq) , 

i> in {z) = Ce lkz u E (k) + De-' kz u E (-k) , (3) 

ipiv{z) = Fe~ qz u E -v {iq) + Ge qz u E - Vo (-iq) , 

^ v {z) =Te lkz u E {k), 



where 



u E (k) 
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(4) 



and 



k = ^^E 2 -m 2 c i 

he 



q = ^m? c ± - {E - Vb) 2 

he 



(5) 

The above coefficients can be determined, as usual, from 
the boundary conditions requiring the wave function to 
be continuous at the potential discontinuities. After some 
simple but tedious algebra we obtain the transmission 
and the reflection amplitudes: 
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cosh(ga) + i — sinh(ga) 
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sinh (qa)e 
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(6) 
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and It is convenient to express the transmission and reflection 

coefficients in terms of their phases as 



R = e i[fe(2a+0-f] ( 1 + a ) 



sinh(ga) 
1 



x { cos( kl) cosh(ga) H (I — a 2 ) sin(fcZ) sinh(qa) 

|_ 2a 

x T, (7) 

where we have introduced R(E) — \R\e itfr (10) 

k (E - V + mc 2 ) . 

a = ; — — . (8) where ip r = (p t — ir /2, while ip t is given by 

q (E + mc/j 

I 



_ 4a(l-a 2 )sinh(2qa)-(l + a 2 ) 2 sin(2fc0[l-cosh(2qa)] ( 

_tan \ 4a 2 [1 + cosh(2ga)] + [1 - cosh(2 ? a)] [(1 - a 2 ) 2 - (1 + a 2 ) 2 cos(2W)] ' ' ( "'' 



Now, accordingly to the stationary phase method, the 
(extrapolated) transmitted and reflected phase times are 
given, respectively, as [l| 



= h 



1 dE 

d(p r 
dE 



E 



E 



- dipt 
l dE 



(12) 
(13) 



Eq 



where we used in these expressions the same central en- 
ergy Eq of the incident wave packet, what is justifiable by 
our previous assumptions about the sharp concentration 
of the initial wave packet around this energy (this corre- 
sponds to the situation in which there is essentially no 
distortion nor reshaping of the transmitted wave, a con- 
dition claimed by several authors as necessary to allow 
a physical meaning to the group velocity 0, [III). The 
above expressions imply the equality of the transmitted 
and the reflected phase times, as it is always the case for 
symmetric potentials [26j. From now on we will refer to 
both these times simply as the phase time t p . Such a 
time corresponds to the (extrapolated) instant in which 
the transmitted and reflected wave packet peak appear 
at z = 2a+l and 2 = 0, respectively. 

Now, b y u sing a general relation obtained by Winful 
et al. [181 127] . we can determine the dwell time r^, which 
is a measure of the time spent by the particle in the po- 
tential region, without distinction of whether it is finally 
reflected or transmitted [H, [2^]. Such a relation, for 
symmetric potentials, reads 

Td=T p -Ti, (14) 

where Tj is the self-interference delay, given by 

(15) 



n 



TYl 



The explicit expressions obtained for the phase and 
dwell time from the above definitions are not particularly 



illuminating and are presented in the Appendix. Here 
we will discuss their properties. The limit of one bar- 
rier (of width 2a) is easily obtained by assuming I = 
and it agrees with the results of [l5[ and Also the 
nonrelativistic limit, obtained by making mc 2 — > oo and 
Vb <C mc 2 , agrees with the results obtained by [||. In 
fact, we can verify these limits directly in the expressions 
for the amplitudes and the transmission phase above. 



III. DISCUSSION AND CONCLUDING 
REMARKS 

Of special interest for us is the limit of opaque barriers, 
qa 3> 1 , in which the phase and dwell times become 



2h 



1 + a 2 
2a 



da\ 
dEj 
m 



2a (k 2 + q 2 ) m 
1 + a 2 k 2 hq 2 



Td = 



(1 + a 2 ) Hq 2 ' 



(16) 
(17) 



where we have used the result that in this limit = 
2a/(l + a 2 )m/(hk 2 ). From the above expressions it is 
clear that both the phase and dwell times saturate in 
the opaque limit, not depending either on the width of 
the barriers or on the distance I between them. This 
demonstrates that the generalized Hartmann effect also 
emerges in the context of relativistic quantum mechan- 
ics. As a consequence, if we extrapolate the concept of the 
group velocity into the potential region, it will be given 
by v g = (2a + l)/r p . This velocity can be made arbitrar- 
ily large, allowing for superluminal group velocities for 
sufficiently large barrier widths a, as we shall see below. 

Figure [5J^ shows the typical behavior of the phase and 
dwell times in the domain of fully relativistic energies, 
where we fixed all the relevant parameters, except the 
barriers width a. We observe that for small values of 
the width a, both the phase and dwell times are greater 
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than the free and the light times. As the barrier become 
thicker, both these times grow slower (in fact, they can 
even decrease, depending on the value of the other pa- 
rameters, as it is the case shown in the figure) and they 
can become smaller than the free and the light times. 
What is more interesting, the phase and dwell times can 
become smaller than the light time even before the satu- 
rated regime is obtained. So, in the domain of fully rel- 
ativistic energies the group velocity can be superluminal 
even before arriving at the opaque limit. A similar behav- 
ior can be observed also in the one barrier case (obtained 
by taking 1 = 0). For comparison, we also showed the 
behavior of the phase time calculated from Schrodinger 
equation. We sec that the nonrelativistic theory predicts 
a phase time of the same order as the relativistic one, and 
sharing the same behavior, specially in what concerns 
the possibility of emerging superluminal group velocities 
before the saturation. Figure 03 shows the same plots 
for the energies in the relativistic scales, but now with 
a greater difference between the energy of the incident 
packet and the height of the barrier. Figures 0\ and 
03 show that Dirac theory can predict group velocities 
which can be smaller (Figure 2A) or greater (Figure 2B) 
than those predicted by the Schrodinger theory, being 
each of these situations determined by the specific choice 
of values for the parameters, especially for the energies. 
These results are in agreement with those observed for 
the single barrier case in reference [lfij], and explain the 
ori gin of the apparent contradictory claims of Krekora et 
al. [151 ] and Leavens and Aers 21 1, concerning to whether 
the relativistic theory predicts group velocities smaller or 
greater than those predicted by the nonrelativistic one. 
In Figure 03 we can check the complete agreement of 
the predictions from both the relativistic and the non- 
relativistic theories in the scale of low (nonrelativistic) 
energies. 

In Figure [3J we plot the behavior of t p and Td for fixed 
a and varying I. Figure [5^ shows, as expected, a trend 
to linear increase with I (except for several resonance 
peaks) as long as a is not very large. That is, outside 
the opaque domain both the dwell and phase times do 
not saturate with the barrier separation I. The same be- 
havior is observed for the nonrelativistic time with the 
same values of the parameters (not shown in the figure) . 
We can also observe the equality between the dwell and 
phase times at resonance (R = 0), as predicted by the 
relation (fl"4")) . Again we observe that the phase time off- 
resonance can be smaller than the light time, even before 
attaining the saturated regime, which imply superlumi- 
nal group velocities. As the barrier width increases, the 
off-resonant phase and dwell times tend to saturate to 
their values at the opaque limit, but still presenting the 
resonant peaks, as we can observe in Figure [3)3. Finally, 
it is only when a — > oo that both these times saturate 
in such a way that the resonant peaks are no longer ob- 
served - the generalized Hartmann effect. Therefore, the 
results of the relativistic theory reinforce the conclusion 
by Winful [j| that the generalized Hartmann effect is just 
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FIG. 2: (Color online) The phase and dwell times vs bar- 
rier width a. We used the system of natural units, in which 
h = c = 1. The energies were expressed in units of the parti- 
cle rest energy m (^ 0.5 MeV for an electron, for example). 
Accordingly, distances and times are given in units of mT 1 . 
For reference we also show the free time in which the peak of 
a free wavepacket traverses the region < z < 2a + I and the 
light time in which a free light pulse traverses the same re- 
gion. For comparison we plotted also the nonrelativistic (NR) 
phase time obtained from the Schrodinger equation [j|. The 
horizontal dashed lines indicate the corresponding saturated 
times in the opaque limit. The first two plots correspond to 
energies in the relativistic scale, while the last one concerns 
to nonrelativistic energies. In all these plots we take I — 0.7. 
(A) E = 1.8, V = 1.5; (B) E = 1.46, V = 2.19; (C) 
Eo = 1.01, Vo = 0.018 (in this last plot the relativistic and 
the nonrelativistic phase times coincide). 



an artifact resulting from taking the opaque limit before 
exploring the variation with I (for an alternative argu- 
ment, using the multiple peaks approach, see [30|). 

Let us now look more carefully into the situation char- 
acterizing the generalized Hartman effect by evaluating 
the flux of particles J z = cip^a z ip both in the region be- 
tween the barriers (as given by ipm), and to the right 
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of the potential for the given energy Eq. The solution of 
the stationary problem gives us 



C 



cosh(qa) + i 



2a 



sinh(<7a) 



Aha 



T 



D = - 



(1 



2a 



3inh( (? a)e lfc(3a+2/) T ! 



for the coefficients in ipm, with T being the transmis- 
sion coefficient. From (j6|) it is plain that in the opaque 
limit T ~ e~ 2qa , so that both C and D decay with bar- 
rier width as e~ qa . Therefore, we conclude that in the 
opaque limit, qa 1, there is essentially no flux, hence, 
no propagation of particles in regions /// and V. Ac- 
cordingly, it follows that the saturated times would be 
the same even if the second barrier were absent [Til ] , sim- 
ilarly to what happens in the nonrelativistic theory [ic| |: 
in fact, expression (jTSJ) is identical to that obtained for 
the relativistic case of a single barrier at the opaque limit 
(see [HI and [HI]). Thus, the condition to the generalized 
Hartman effect to occur is the condition of no transmis- 
sion, in which case it makes no sense to associate any 
velocity to the tunneling process (ij. 

On the other hand, it is possible to have situations in 
which qa is large, but finite, such that there is still an 
appreciable transmission (before the saturation regime), 
and the associated group velocities during the tunnel- 
ing are supcrluminal. Notice that we have considered 
wavepackets sharply centered around a given energy in 
such a way that the transmitted wave packet could be 
seen essentially as a (attenuated) nondistorted version of 
the incident one, a feature that is claimed by several au- 
thors as allowing one to attribute a physical meaning to 
the group velocity [25[ Therefore, if one maintains 
the interpretation that the group velocities are propa- 
gation velocities, it would seem that relativity does not 
forbid supcrluminal tunneling velocities in the single or 
double barrier tunneling. However, it must be noticed 
that while there is little doubt that the group velocity in 
the region V (or region I, for that matter) has a physical 
meaning, superluminal group velocities emerge when we 
extrapolate the concept of group velocity to the region 
within the barriers. But, it is clear that inside the barri- 
ers the (evanescent) wavepacket undergoes great distor- 
tion, not sharing the same shape as the incident or trans- 
mitted ones; in fact, within the barriers the wavepacket 
does not even have a pea k which travels from one bound- 
ary to the another [ilJlClll5j|. Thus, despite the fact that 
the group velocity has well defined meaning for the in- 
cident (before reflections) and transmitted regions, the 
extrapolation of this concept to the region inside the 
barriers cannot be justified, and consequently there is 
no justification for associating this (extrapolated) group 
velocity with tunneling velocity. 

In what concerns the dwell time, since it does does not 
distinguish between the transmitted and reflected chan- 
nels, it is better interpreted as a cavity lifetime. However, 
it is important to notice that here, contrary to what hap- 
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FIG. 3: (Color online) The phase and dwell times vs sepa- 
ration I between the barriers. Eo ~ 1.8 and Vo = 1.5. (A) 
a — 0.7. (B) a = 3.0. The horizontal dashed lines indicate 
the corresponding saturated times in the opaque limit. We 
used the same (natural) system of units as in the Figure 2. 
The peaks in both figures correspond to resonant tunneling. 



pens for a Fabry-Perot cavity 0, [T(| , the phase and dwell 
times are not equal in the off-resonance case [see Figs. 2 
and 3 and the limits (fTo| - (|17[) ]. This prevents an imme- 
diate identification of phase time as a cavity lifetime in 
the present scenario. 

Summarizing, in this paper we have analyzed the rela- 
tivistic tunneling of a spin 1/2 particle through two suc- 
cessive electrostatic potential barriers and showed that 
the so called generalized Hartmann effect also occurs in 
the realm of relativistic quantum mechanics. In addition, 
we obtained that the dwell time also saturates with the 
width of the barriers. We demonstrated that the phase 
and dwell times can become smaller than the light time 
(which implies superluminal group velocities) even be- 
fore the saturated regime is obtained, and we observed 
that the group velocities predicted by the relativistic the- 
ory can be smaller or greater than those predicted by 
Schrodinger's theory, depending on the values of the pa- 
rameters. We also showed that the phase and dwell times 
show an almost linear increasing with the separation be- 
tween the barriers, and tends to saturate only when the 
barrier becomes extremely opaque. Finally, we discussed 
a possible interpretation of the results, favoring the ar- 
gument that the group velocity cannot be interpreted as 
a tunneling velocity. 
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APPENDIX 



1 



(kl) 



E I h x 

fe2 ~ k 2 q 2 T 2 + A 2 



In this Appendix we list the explicit expressions for the 
phase time, r p , and the self-interference delay, Tj. From where we have defined 



(A.l) 



and 



T = 8a 2 cosh(2ga) - 4(1 + a 2 ) 2 sin 2 (kl) sinh 2 (ga) , 
A = 4a(l - a 2 ) sinh(2qa) + 2(1 + a 2 ) 2 sin(2kl) sinh 2 ((?a) , 

fti = A {2(1 + a 2 ) [(1 + a 2 )Eq 2 (2kl) sin(2fcZ) - 4a 2 mc 2 (fc 2 + q 2 ) cos(2fc/)] sinh 2 (qa) 

- 4a 2 mc 2 (fc 2 + q 2 ) [(1 + a 2 ) + (3 - a 2 ) cosh(2ga)] 

+ k 2 (2qa)(E - V ) [(1 + a 2 ) 2 cos(2A;Z) - (1 - 6a 2 + a 4 )] sinh(2ga)} 

+ T {-4a(l - a 2 )k 2 (2qa)(E - V ) cosh(2ga) 

+ 2(1 + a 2 ) [(1 + a 2 )Eq 2 (2kl) cos(2W) + 4a 2 mc 2 (fc 2 + q 2 ) sin(2jfe/)] sinh 2 (qa) 

+ [4a(l - 3a 2 )mc 2 (fc 2 + q 2 ) - (1 + a 2 ) 2 k 2 (2qa)(E - V ) sin(2jfe/)] sinh(2ga)} . 
From equations ((6]) and (|7|), the self-interference delay, as defined in (|15| . is given by 

m (1 + a 2 ) hi 



hk 2 4a 3 h 3 



with 



1 



hi = 2 a (^ — of) sin(2fc£) sinh (2qa) + a cos (kl) sinh(4ga) 

+ a(l - a 2 ) sm(2kl) sinh 2 (ga) cosh(2ga) + (1 - a 2 ) 2 sin 2 (£;Z) sinh 2 (^a) sinh(2ga) , 
1 



8a 4 



{8a 4 cosh 4 (e?a) + [l + 6a 4 + a 8 - (1 - a 4 ) 2 cos(2fcZ)] sinh 4 (qa) 



(A.2) 
(A.3) 



(A.4) 
(A.5) 



(A.6) 



+ a 2 [(1 - a 2 ) 2 + (1 + a 2 ) 2 cos(2W)] sinh 2 (2ga) + 2a(l - a 2 )(l + a 2 ) 2 sin(2fe0 sinh 2 (ga) sinh(2ga)} . (A.7) 

Finally, the dwell time is obtained from the phase time and the self-interference delay from (fl"4"|) . 
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